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, We present an extension of the modular recursive Green's function method (MRGM) for bal- 

' listic quantum transport to include magnetic fields. Dividing the non-separable two-dimensional 

scattering problem into separable substructures allows us to calculate transport coefficients and 
, scattering wavefunctions very efficiently. Previously unattainable energy and magnetic field regions 

can thereby be covered with high accuracy. The method is applied to magnetotransport through a 
circle and a stadium shaped quantum dot at strong magnetic fields and high energies. In the edge 
state regime we observe strong multi-frequency Aharonov-Bohm oscillations. By analyzing them in 
terms of a multi-channel interference model, we classify these fiuctuations within the framework of 
Fano resonances and discuss their geometry independence. For high energies (mode numbers) we 
observe localization of the scattering wavefunction near clcissical trajectories. 
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I. INTRODUCTION 



Accurate simulations of ballistic transport through quantum dots have remained a computational challenge, despite 
• the conceptional simplicity of the problem. This is in part due to the fact that many of the most interesting phenomena 
' occur in a parameter regime of either high magnetic field B or high Fermi energy Ep- The regime of strong magnetic 
I field B, where the magnetic length (in a.u.) Ib = \J cjB is small compared to the linear dimension D of the dot, 
■ O Ib D, gives rise to the Quantum Hall effect the Hofstadter butterflyj^ and Aharonov-Bohm oscillations of transport 
^ coefficients f2. The high energy domain, where the De Broglie wavelength Xd = V^Ep satisfies Ad ^ D, is of particular 
relevance for approaching the semiclassical limit of quantum transport and for investigations of quantum signatures of 
classical chaos^i^ Both of these regimes pose considerable difficulties for a numerical treatment. Methods based on 
the expansion of the scattering wavefunction in plane or spherical waves become invalid at high fields since diamagnetic 
contributions are generally neglectedi^ Methods employing a discretization on a grid are limited by the constraint 
that the magnetic flux per unit cell must be small, which, in turn, requires high grid densities for large B.^ The same 
' requirement has to be met for high Ep, where many grid points are needed to accurately describe the continuum limit. 
] This implies however a large number of inversions of high-dimensional matrices and therefore limits the applicability 
of these approaches for large B and (or) large Ep. 

In the current communication we propose an approach that allows accurate treatment of these regimes. We present an 
, extension of the previouslyS, introduced Modular Recursive Green's function Method (MRGM) to include an additional 
magnetic field perpendicular to the two-dimensional scattering surface. The underlying idea for our approach goes 
back to Sols et alM and to the widely used Recursive Green's Function Method (RGM)i2iii In the standard RGM 
the Green's function is propagated through the scattering region from one transverse strip to the next by repeated 
H ' solutions of a matrix Dyson equation. We show that the efficiency of this conventional discretization can be increased 
considerably by taking the symmetries of a scattering problem into account from the outset. Specifically, when the 
two-dimensional nonseparable open quantum dot can be built up out of simpler separable substructures (referred to 
O 1 in the following as modules), the Green's functions for each of these modules can be calculated efficiently and virtually 
J-J ' exactly. Calculation of the S matrix and of the scattering wavefunction is thus reduced to "welding" together the 
^ , modules by a very small number of recursions. Key to this approach are tight-binding grids which are symmetry- 
adapted for each module. This leads to the separability of the eigenf unctions in the modules and allows an efficient 
incorporation of boundary conditions. As a result, a much higher grid density can be easily handled, which, in turn, 
' is prerequisite for treating short magnetic lengths Ib and short wave lengths Ad . Matrix Dyson equations have to be 
■ ■ ' solved only for each junction between the modules. The total number of necessary recursions (i.e. high-dimensional 
matrix inversions) is thereby reduced to the number of modules needed to reconstruct the quantum dot. 
The efficiency of the MRGM will be demonstrated by applying it to transport through a circular and a stadium 
shaped quantum dot. These systems are known as prototype structures for regular and chaotic dynamics and have 
been studied thoroughly in the literature A^^Sii^ Concerning the theoretical approaches for the investigation of electron 
dynamics in quantum dots, considerable attention has been dedicated to reach higher energie9i2*i^iiS*i& and higher 
magnetic fieldsiiSiiLiSiiSiSS Especially for the study of transport through open stadium billiards, several different 
methods have been employed iiSiSLSSiSii'* In the following we will present numerical results obtained by the MRGM 
which attain a parameter range, to our knowledge not yet explored by other approaches. For small Ad we investigate 
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the localization of the scattering wavefunction near classical scattering trajectories. Characteristic differences in the 
dynamics of generically regular and chaotic systems will be highlighted. In the high magnetic field regime, which 
is governed by edge states, differences between the dynamics in different geometries disappear and are replaced by 
universal quasi-periodic conductance oscillations. At a critical magnetic field these oscillations break off and transport 
terminates entirely. In the regime where more than one edge state is excited in the dot, we find interference fluctuations 
which we analyze in terms of a multi-channel Fano interference model. The key to the understanding of the observed 
fluctuations is that inter-channel scattering between different edge states takes place only by diffractive scattering at 
the lead junctions. 

This paper is organized as follows. In section II the method for inclusion of a magnetic field in the MRGM is 
presented. Section III is dedicated to a discussion of numerical results, illustrating the high magnetic field and high 
energy behaviour in quantum dots. The paper concludes with a short summary in section IV. 



II. METHOD 



We consider ballistic nanostructures with a constant electrostatic potential inside the two-dimensional cavity, impose 
hard- wall boundary conditions, and assume a constant magnetic field to be oriented perpendicular to the scattering 
plane. The shape of the quantum dot will be chosen to be either a circle or a stadium (see Figs. I5I6I below), which 
represent prototype systems for regular and chaotic classical dynamics, respectively. Two semi-infinite waveguides 
of width d at different electrochemical potentials (/ii, /Z2) are attached. The aperture of the leads is chosen to be 
very smaU d/D = d/^/A^ = 0.0935, where A'^°^ = 4 + TT is the scaled area of all the cavities studied and _D is a 
characteristic linear dimension of the cavity. At asymptotic distances, i.e. far away from the quantum dot, scattering 
boundary conditions are imposed. The asymptotic scattering state can be factorized into a longitudinal flux-carrying 
plane wave and a transverse standing wave. The latter is a simple sine wave in the fleld-free case and a combination 
of Kummer functions when the magnetic field is turned on i^^i^^ In our local coordinate system the longitudinal 
(transverse) direction in the i-th lead is always denoted by Xi (yi). The wavefunctions in the waveguides thus vanish 
at Hi = ±d/2. Atomic units {H = \e\ = rrics = I) will be used from now on, unless explicitly stated otherwise. 



A. Brief review of the MRGM 



In order to highlight the technical difficulties in incorporating a magnetic field we start by briefiy reviewing the 
MRGM for the field-free case. Starting point is the observation that a large class of dot geometries with non-separable 
boundaries can be decomposed into separable two-dimensional substructures, referred to in the following as modules. 
For each of these modules the discretization of the corresponding tight-binding (tb) Hamiltonian can be performed 
on a symmetry-adapted grid. The grid for each module is chosen such that the eigenfunctions of the tb Hamiltonian 



77"^ = ^e,|0(*l+E^^^I*)0'l (2.1) 
separate into two generalized coordinates. H^^ contains hopping potentials Vij for nearest-neighbour coupling and site 
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energies Ei. Both quantities are chosen such that the Schrodinger equation, H^^\ipm) = -E'mlV'm), converges towards 
the continuum Schrodinger equation in the limit of high grid point density. The most straightforward application 
of this approach refers to modules for which the boundaries are nodal lines of Cartesian {x,y) or polar coordinates 
{g^f). For these coordinate systems we haveS at i? = 



rre _ -ei±l/2 TrV _ -1 ^ _ _A \ L 



i,i±l — 2QiAQ^ ' " i.j±\ ~ 2gfAip^ ' ~ Ag^ T g'fAip'^ 



(2.2) 



with Qi — \i — l/2\Ag. For separable energy eigenfunctions of the general form — \Ek) <E) \Ek.n) the spectral 

representation of the retarded (-I-) and advanced (— ) Green's function G*(r, r', B, Ep) of the module is simply given 
by 

G^(r,r-,i3,i^,) = E(a|i^.)(i^.|c.0E T+"^^'''"/^^ - (^.3) 

k n 

where (a, (3) stand for the (generalized) coordinates {x, y) or ((/>, p). The indices (fc, n) represent the quantum numbers 
of the separable eigenfunctions \Ek)^ \Ek^m) associated with the degrees of freedom a and (3 respectively. 
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The Green's functions of the separate modules are joined by solving a matrix Dyson equation, 

G = G" + G"VG, (2.4) 

where G'^ and G denote Green's functions of the disconnected and the connected modules, respectively. The matrix 
V denotes the hopping potential V multiplied by the size of the unit cell V = VAr, which in a Cartesian (polar) grid 
is Aji = Ax Ay {= giAgAip). The complete scattering structure can thus be assembled from the individual modules 
(much like a jigsaw puzzle). The number of necessary recursions [i.e. solutions of (|2.4|) ] is (approximately) equal to the 
number of modules. The exact number depends on the number of link modules required for different grid structures. 
For example, in order to connect a half-circle with a rectangle we need one additional link module which is plugged in 
between [see Ref. |3 for details] . The key property of these link modules is their adaption to two grid symmetries [see 
Fig. IJh] ■ Mathematically speaking, the transition from a polar to a Cartesian grid requires a link module in order to 
preserve the hermiticity of the tb Hamiltonian at the junction. In the recursion the link module is connected to the 
Cartesian (polar) grid by means of the hopping potential V^{V'^), respectively. 

Once the Green's function G^ for the combination of all modules is assembled, the transmission amplitudes tnm from 
entrance lead mode m into exit lead mode n can be calculated by projecting G"*" onto the transverse wavefunctions 
in the leads XniUi)- With the corresponding longitudinal wave numbers kx^.n we have (at zero magnetic field), 

/d/2 i-d/2 
dy2 dyi Xniy2)G+{y2,yi,EF)xrniyi) ■ (2.5) 
-d/2 J-d/2 

Together with the reflection amplitudes r„m (for which an analogous relation holds) the S'-matrix is completely 
determined and satisfies the unitarity condition implied by current conservation, 

M 

^(|tnm|' + |r„™|') = l. (2.6) 
n=l 

The integer M denotes the number of open channels in the leads. According to the Landauer formula, the total 
conductance g through the quantum dot is given by 

1 1 

5 = - V |i„„|2 = - r'°* with T*°' + = M . (2.7) 

m.n—l 



B. Inclusion of the magnetic field 



Incorporation of the magnetic field into the MRGM poses a number of complications. The solutions of these 
difficulties will be presented in this section. At the core of the problem is the preservation of separability of the 
Schrodinger equation. The usage of gauge transformations as well as of Dyson equations for decomposing non- 
separable structures into separable substructures plays a key role in accomplishing this goal. The field B = (0, 0, B) 
enters the tb Hamiltonian H2.1II by means of a Peierls phase factor 



Vr. 



r,r' • exp 



i/c / A(r)dr 



(2.8) 



with which the field-free hopping potential Vry is multiplied. The vector potential A(r) satisfies V x A(r) = B. The 
Peierls phase will, of course, in most cases destroy the separability of the eigenfunctions of H^^. The difficulties can 
be, in part, circumvented by exploiting the gauge freedom of the vector potential, i.e.. 



A' 



VA, 



(2.9) 



where A(r) is a scalar function. By an appropriate choice of A the wavefunction remains separable on a given symmetry 
adapted grid. Specifically, to preserve separability we employ the Landau gauge for a Cartesian grid 



A = (-By, 0,0) , 
and the "symmetric" or circular gauge for a polar grid 

A = B/2{-y, x,0)^ px B/2 . 



(2.10a) 
(2.10b) 
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The scalar gauge potential generating the gauge transformation from H2.10a() to (|2.10b|) is X{x, y) = Bxy/2. 
A major complication results from the fact, that in the presence of the magnetic field the separability on an unrestricted 
grid of a given symmetry does not imply the separability in the presence of boundary conditions of the same symmetry. 
We illustrate this problem with the help of one typical example, the semi-infinite quantum wire with lead width d 
(Fig. We impose hard-wall boundary conditions ip{x,y — ztd/2) — and consider first the infinite quantum wire 
along the x direction. Because of the Cartesian boundary conditions, the symmetry adapted gauge is the Landau 
gauge A = —Byk. Consider, for notational simplicity, the Schrodinger equation in the continuum limit, 

H<^iy. .) 4 (p + i a) ^ Hx. y) = l{-^-$- ^4 + ^) ^) - ^'-''^ 

Since the longitudinal momentum px = --id/dx commutes with H , the separability of the wavcfunction persists in 
the presence of the magnetic field: (t>{x^y) = fk{x)x{y) with fk{x) — e**^^. If, however, one introduces an additional 
Cartesian boundary condition along the y-axis [i.e. = 0,y) = for a semi-infinite lead] the situation changes. In 
the absence of the magnetic field, B = 0, the linear term in vanishes and thus the choice f{x) = sin(fca;) [i.e. a 
linear combination of /±fc(a;)] satisfies the boundary condition and preserves the separability, even though p^ is no 
longer conserved in the semi-infinite lead. However, for B ^ and the same boundary condition ^{x — 0,y) = 0, the 
term linear in B and p^ destroys the separability. The wavcfunction takes now the general form 

(^(x, y)^Yl c^nXmn{y). (2.12) 

m n 

The breakdown of separability by the introduction of an additional boundary condition indicates that the Green's 
function of confined modules will be more complex than for extended systems for the same symmetry adapted grid 
and the same gauge. Therefore, the program of the modular method of building-up extended complex structures by 
"welding together" smaller modules of higher symmetry will be executed in reverse: non-separable confined modules 
will be generated by "cutting in pieces" larger separable modules. Confining boundary conditions will be introduced 
rather than removed by the matrix Dyson equation. In the example above, the semi-infinite quantum wire is generated 
by cutting the infinite wire at the line a; = 0, thereby imposing the additional boundary condition. Just as connecting 
modules, so is disconnecting a given module equivalent to the application of a matrix Dyson equation, 

G'^ ^G^ + G^VG^. (2.13) 

In this context G^ {G^) is the Green's function of the extended (confined) module and V is the hopping potential 
that connects the modules. Solving (|2.13l) in reversed mode (i.e. for G^ rather than for G^) amounts to dissecting 
the larger module. 

Provided that the Green's functions of all the necessary modules are available, we have to link them with each other 
to assemble the entire scattering geometry. However, in the presence of a magnetic field we have to take into account 
that the different modules will be calculated in different symmetry-adapted gauges. Joining different modules requires, 
therefore, in general a gauge transformation. For the Green's function on the grid G'(ri,r^) this transformation is 
simplified by the fact that the matrix of gauge transformations 

[A(r,)]^., -exp[-iA(r,)/c]J,fe (2.14) 

is diagonal in the grid representation. Correspondingly, the transformation of both the hopping potential V and the 
Green's function is local, i.e. 

F(r„ r;.) ^ \/'(r,, r^) = A(r,) V{v,,v'^W{v'^) (2.15) 



G(r„r;-) G'(r„r;-) = A(r,) G(r„ r;-)A*(r;-) . 

It is thus not necessary to transform gauges of different modules to one global gauge. Instead, it is sufficient to 
perform a local gauge transformation at the points of the junctions {r^}, such that the gauges of the two modules to 
be joined agree at these points. 

Finally, in order to extract the 5- matrix, i.e. the amplitudes t„„i and r„m, matrix elements of the current operator 
must be of gauge-invariant form. This requirement can be fulfilled by employing a double-sided gradient operator 
which is defined aaSS 



ffig = /(x)D5(x) - g(x)D7(x) = -gf)f with D = V - -A(x) 



(2.16) 
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With its help the transmission amplitudes can be evaluated asiSiSSiSS 

'iy^n^m J-d/2 J-d/2 

G+(x2,x;,i?j.,S)(D'.X'i)x,„(2/'i)e^''"^'i . (2.17) 

The unit vectors x„ are assumed to be pointing in outward direction of the n-th lead and 6m denotes the outgoing 
particle flux carried by Xm(yi)e*'^"*^i through the lead cross section. Determination of transverse states XmiUi) and 
of the corresponding longitudinal momentum km as well as the normalization factors 9m will be discussed below. For 
reflection amplitudes Vnm, a relation similar to H2.17|l holds. From tnm and r„m the conductance can be calculated 
by means of the Landauer formula [Eq. (|2.7|) ]. 



C. Calculation of modules 



This section is dedicated to the evaluation of the Green's functions for those modules, which we need to assemble a 
circle and a stadium billiard: the half-infinite leads, the rectangle, the circle, and the half-circle. With the exception of 
the circle, for all these modules Eq. (|2.3|) is not applicable. This is due to the non-separability for confined geometries 
as discussed above. Moreover the spectrum in open structures like the semi-infinite lead is continuous rather than 
discrete. Unlike in the field-free case^^ the resulting integrals cannot be calculated analytically. However, both 
problems can be overcome by applying the matrix Dyson equation in a non-standard way. 



1. Rectangular module 

As illustrated above for the semi-inflnite wave guide, the Dirichlet boundary condition for the confined structure 
of a rectangle with magnetic field is not separable, no matter which gauge is chosen. The separability can however 
be restored by imposing periodic boundary conditions on two opposing sides of the rectangle. Topologically, this 
corresponds to folding the rectangle to the surface of a cylinder (Fig. In this case we connect the flrst (P) and 
the last (Q) transverse grid shce of a rectangular grid by a hopping potential |V^q| = \Vqp\ = -rf^- The Green's 
function of this "cylinder surface" (cs) will be denoted by C^^ in the following. The calculation of the rectangle Green's 
function G" will be obtained out of G"^^ by a Dyson equation used here in "reversed" mode, i.e. for disconnecting tb 
grids. This method for calculating the rectangular module may seem like a detour, but it is numerically more efficient 
than a strip-by-strip recursion. For completeness we mention that an alternative way to calculate was proposed 
in Ref.m 

The Green's function for the cylinder surface C^^ can be constructed from separable eigenfunctions, \Em) = l-E^) ® 
l-^fcn)i according to Eq. 1)2.3(1 . Solving the tight-binding Schrodinger equation for the cylinder surface, we obtain 
for the longitudinal eigenstates = (A^^r Ax)^^/^ exp(i27rfcj/7V2;), which results in a tridiagonal, symmetric 

matrix-eigenproblem of size Ny x A^j^ for the transverse modes^^ 



Ekn{yi\El^) - ^ 



,2'Kk B ^ , 



>^{yi\EL) - ^{{m-i\EL) - ^imlEL) + im+ilEL)). (2.18) 

By "cutting the cylinder surface open" along a line of constant x, we obtain from G'^'^ the desired Green's function 
G"^ for the rectangle (Fig.|3i). We demonstrate this for the rectangle Green's function Gp^ from the first transverse 
slice P to any other slice X. To determine Gp^^ we solve the following system of Dyson equations, 

Gpx = Gfx - GpqVQpGfx - Gpp^pqGqx (2.19) 
GpQ = C!pQ + GpqVqpGpq + GppVpqGqq (2.20) 
Gpp = Gpp + G^pqVqpGpp + G^ppVpqGqp , (2-21) 

where the first line is the "reversed" Dyson equation. The three unknowns in the above equations are the Green's 
functions connecting the slices {P,X),{P,Q) and {P,P), Gpx,GpQ,Gpp. By solving these three equations, the 
unknowns can be uniquely determined. 
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2. Circle and half-circle 



In symmetric gauge, A = B /2(^y, x,0), the Dirichlet boundary value problem for the circle with magnetic field 



is separable, \E„ 



\E 



kni 



On a discrete tb lattice this statement remains true, provided that a circular 



grid is employed. With the eigenstates for the azimuthal degree of freedom, 
and radial eigenstates gkn{Qi) = ^/ol x (ftl^-fe 
symmetric eigenproblem 



{^,\El) = (7V^A(^)-i/2exp(i2^A:j77V^) 
E^„), the finite difference equation for the gkniffi) results in a tridiagonal 



Ekn gkniSi) 



1 



1 

2A^ 



cos I — -S^— — - - 1 

2c / 



gkniQi) 



■.gkniQi^l) - 2gkn{Qi) H ^ , gkn{Qi+l, 



(2.22) 



The Green's function for the circular module is then calculated by a straight-forward application of Eq. (2.3). 
For the Green's function of the half-circle we employ an analogous procedure as in the previous subsection: we dissect 
the circle Green's function into half-circles by means of a "reversed" Dyson equation. We demonstrate this by way of 
the example depicted in Fig. [St, where the "full circle" (fc) is split up into two "half-circles" (he). The resulting two 
halves are almost identical, with the exception of the two additional radial grid slices, by which the right half-circle 
is larger. For assembling the stadium billiard we have to make sure that the tb grid of the half-circle module can be 
linked directly to a vertical grid [see Fig.l^b]. For this reason, only the left one of the two half-circles in Fig.[5lD can 
be used for this purpose. 

Consider as example the Green's function G"^^ describing the propagation from the grid slice P at the junction of 
the two half-circles to any radial grid slice X situated on the "left half-circle" (see Fig. |2b) ■ G^x is determined by 
the following system of Dyson equations 



Ghc y^fc y^hc T/ y-^fc 

— ^PX ~ '-'PP^PQ^QX 



G 



PX 
fc 

PP 



Ghc I /-l\\C f^fc 
PP -\- Upp ypQ^Qp 



(2.23) 
(2.24) 



which yields a unique solution for G^x ■ 



3. Semi-infinite lead 

Because of its continuous spectrum, the Green's function for the semi-infinite lead poses an additional challenge 
beyond that of the non-separability of the wavefunction discussed above. We therefore apply one further "trick" to 
bypass this problem. Our approach is based on the observation that adding a slice to a semi-infinite quantum wire 
leaves this wire (up to irrelevant phases) invariant (see Fig. We assume a semi-infinite lead with x € [Ax,oo) 
and hard-wall boundary conditions at x = Ax and y = ±(i/2. To this object we add a slice consisting of just one 
transverse chain of tb grid points which we place at a; = 0. The system of Green's functions for the propagation from 
the transverse chain at x = (P) back to itself {P) or to the first transverse slice of the semi-infinite lead [Q) at Ax 
reads 

Gpp = Gpp -\- GppVpgGgp, (2.25) 
Gqp — GqqVqpGpp. (2.26) 

Each multiplication involves a matrix product with a dimension equal to the number of transverse grid points. The 
key point is now that the system of Eqs. (|2.25l2.2(j|) can be closed through the invariance condition (Fig. |21i) for the 
semi-infinite lead, i.e. Gpp = Gqq. In Landau gauge A = {—By, 0,0) the latter relation does not involve additional 
gauge phases since these are already contained in the hopping matrix element. We additionally note that an equivalent 
point of departure for the derivation of Gpp is the Bloch condition for states in the lea,d}^'^^ 

Setting Z = GppVqp and using the hermiticity condition Vqp ~ Vpq = V* , Eqs. H2.25l2.2t)|l can be converted to a 
quadratic matrix equation 

ZZ -V-'^{G°pp)-^Z + V-'^V* = 0. (2.27) 

Solvents Z of a quadratic matrix equation Q{Z) = can be constructed from the eigenpairs {(3i,Xi) of the corre- 
sponding quadratic eigenvalue equation Q{l3i)xi — 0, i [1, . . . , 2N] in the diagonal form,3i 



Z^MBM-^ with M= [xi,...,Xiv], ^ = diag(/30- 



(2.28) 
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The quadratic eigenvalue equation is equivalent to a generalized eigenvalue problem Ax = PCx of twice the original 
dimension. '^^ Its 2N dimensional eigenvectors x — (x? Px) solutions of the symmetric eigenproblem 



V J \ PxJ ' \ V J \ fx 

where (Gpp) ^ — Ep — and is the Hamiltionian of the one-dimensional transverse tb strip at x = 0. The 
Fermi energy Ep and the magnetic field B enter (I2.29|) as independent parameters at which the eigenstates Xm and 
eigenvalues (3m are evaluated. The longitudinal momenta of the lead states ^rn{x,y) — Xm{y)e'^^™'^ / V^m are related 
to the eigenvalues by the relation (3 — exp{ikAx). The orthonormalisation and the completeness relations of the 2N 
eigenvectors Xm can be formulated in terms of matrix relations, for the generalized eigenproblem, 

xlCxn^2^-^5mn and y^^^ = 2^-^n-\ (2.30) 



y/dmOn ™ l^ml ^ 

With this specific choice of normalization the norm factors Qm are determined such that every propagating state carries 
unit flux. We note parenthetically that the quadratic eigenvalue equation could also be applied to the semi-infinite 
lead at zero B field. However, in that case, the Green's function for quantum wires can be calculated analyticallp^^ 
by complex contour integration. 



4. Scattering wave functions and efficiency of the MRGM 

The MRGM is particularly well-suited to determine transport coefficients as the Green's function is then required 
only at the junctions between the modules and does not have to be evaluated throughout the interior of the entire 
quantum dot. Since for the calculation of the scattering wavefunction the Green's function throughout the entire 
scattering region is needed, this particular advantage cannot be made use of here. However, also in this case, the 
MRGM is still more efficient than the standard RGM, as will be explained below. 

The wavefunction '(/'(x) can be obtained at any point x by projecting the retarded Green's function (by means of the 
operator D) on the incoming wave (in mode m)^Si2S 

V'™(x) = --^ f'^ dy[G+{^,^\Ep,B) (D' • x'J x™(2/'i)e*"""'^ ■ (2.31) 

^yUm J-d/2 

contains the solution of the Dyson equations for all linked modules. That the evaluation of Eq. (|2.31|) can be done 
very efficiently results from two properties: First the number of recursions (i.e., of matrix inversions) needed to obtain 
G^ is given by the fixed number of modules required to build up the scattering structure. This number is independent 
of the De Broglie wavelength. The latter enters only in terms of the size of the matrices involved in the recursion, 
since with increasing Ep (decreasing Xd) more grid points are required to represent the continuum limit. Compared 
to the standard RGM the numerical effort is reduced since in that approach the number of recursions scales with the 
grid density, i.e. oc kp. A second advantage of the MRGM concerns the incorporation of the boundary conditions. 
In the modular method the boundaries follow by construction the nodal lines of the symmetry-adapted coordinate 
system for the module. Due to this reason the convergence towards the continuum limit is enhanced as compared to 
the slice- by-slice recursion. The calculation of the transport coefficients as a function of the Fermi wavenumber kp (or 
Fermi energy Ep) is simplified by the fact that the solution of the eigenvalue problem {\Em), Em) entering the Green's 
function for each module [Eq. H2.3I) ] is independent of Ep. For the evaluation of the Green's function at different 
values of Ep the eigenproblem W^^\Em) — Em\E„i) therefore has to be solved only once. Unfortunately, this feature 
does not extend to the variation of the magnetic field since both !£",„) and _E,„ are dependent on B. Because of this 
property a new solution of the tb eigenproblem is required for each value of the field. The most severe restriction of 
the MRGM is, however, that its applicability is limited to those scattering structures which can be assembled from 
or cut out of separable modules. Also random potentials and soft walls can only be included as long as they preserve 
the separability of each module. We mention at this point, that a "hybrid RGM" for dealing with arbitrary boundary 
geometries was presented in the literature 



III. NUMERICAL RESULTS FOR HIGH B FIELDS AND LARGE kp 



In this section we present first magnetoconductance results which were calculated within the MRGM at high 
magnetic fields B and large Fermi wavenumbers kp. As prototype shapes of the cavity we use the circle and the 
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Bunimovich stadium and consider different geometries for the attached quantum wires. 

A. Accuracy checks 

Several checks for the accuracy of the numerical results have been performed. Exact relationships for transport 
coefficients such as conservation of unitarity and the Onsager relations are fulfilled with an accuracy of better than 
10^^*^. The grid density is chosen such that the magnetic flux per unit cell is BA^/c < 0.01 (as in Ref . Is^h . Moreover, 
the typical number of grid points per Fermi half- wavelength is greater than 30. Only for very high energy calculations 
(Fig. ISJ the relative grid density is lower. For low magnetic fields, we can compare our results for |tnm(fcF)P with 
previous methods. As an example we show in Fig. 01a comparison for |tii(fci?)P with the calculation by Yang et al.^ 
which is based on a wave function expansion in spherical waves. The agreement for the circle is very good although 
diamagnetic terms are neglected in the approach of Ref. 0. For the stadium, the differences between the two methods 
are somewhat larger. This is due to the fact that the expansion of the stadium wave functions in spherical waves leads 
to a unitarity deficiency (see Fig. 0)). We can also reproduce previous results of Ref. concerning weak- localization 
line shapes and statistical magnetoconductance properties in chaotic and regular cavities. These will not be treated 
again. Our focus will be on the high magnetic field and high energy regime, where other methods failed. 

B. Wavefunctions 

The starting point for the analysis of the scattering states V'(^) fo^' ballistic transport through quantum dots is 
Eq. (|2.31|) . Figures 13 and display the resulting electron density oc |V'(r)P in the scattering region. In Fig. [S] 
we consider the wavefunctions at very high kp for both the circle and the stadium billiard, which are prototypical 
structures for regular and chaotic dynamics respectively. Large kp corresponds to the regime where the convergence 
towards classical scattering trajectories is expected to emerge. Figures and[S|3 illustrate the different dynamics 
for an injection at high (m = 20) and at low mode number (m = 1), respectively. Since high mode numbers 
correspond classically to a large injection angle, the wavefunction condenses around a pentagon-shaped whispering 
gallery trajectory. For low-mode injection, a small circle representing the centrifugal barrier (or caustic) is seen, as 
well as rays representing the asterisk orbits.'^^ Figures l^; and display scattering states for the stadium. At low 
magnetic fields, the dynamics is chaotic and a typical wavefunction features a quasi-random pattern with a modest 
density enhancement near classically unstable periodic orbits (not shown). For special configurations of fc^? and B 
"scars" emerge in the scattering wavefunctions (Fig. |S|;). By contrast, for high magnetic fields the classical motion 
becomes regular. In the present example the wavefunction condenses around a "bundle" of cyclotron orbits executing 
three bounces at the cavity wall before exiting by the entrance lead (Fig.|3i). There has been an extensive discussion 
in the literature as to the existence of scars in open quantum billiards.^^ Our present results clearly underscore that 
scars, defined here as the condensation of the wavefunction near classical (not necessarily unstable) trajectories, clearly 
exist for large kp- Figure illustrates the formation of edge states at high fields. With increasing B fewer edge states 
can be excited in the cavity. In Fig.E): {B = 68.5) three transverse edge states are present while in Fig.|SJl {B = 125) 
only a single edge state remains. For two edge states carrying flux across the quantum dot, interferences give rise to 
a stationary nodal pattern with a fixed number of antinodes along the boundary (see Fig.|B^,b). We note that we are 
not aware of any other method that has so far been capable of investigating scattering states of open structures in 
this high-magnetic field regime. 

C. Transport coefficients 

The interference between different edge states gives rise to characteristic fluctuations in the transport coefficients. 
Figure shows the high magnetic field regime of the transmission probability in the first mode m = n = 1 for both 
circle and stadium. Different orientations of the exiting quantum wire were chosen (oriented 90° and 180° relative to 
the incoming lead). A few universal trends are easily discernible: above a certain critical value of the magnetic field 
(denoted by B]), the strongly fluctuating transmission probability gives way to very regular oscillations in all four 
cases [see insets of Fig. |7|for magnification]. The threshold value B]. and the magnetic field, at which transport is 
terminated (separately displayed in Fig. IS)) are identical for all systems investigated. Below B]. the transport signal 
displays "beats" , i.e. the Fourier transform of the signal is characterized by several frequencies. The "universality" 
(i.e. geometry independence) of these features is related to the fact that in the high magnetic field regime transport is 
controlled by edge states (as depicted in Fig.E)). These states play a very prominent role in the Quantum Hall effect 
and have been studied extensivelyii*ii*2£i2& At magnetic fields, where the magnetic length is smaller than the system 
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dimensions, Ib ^ D, they are the only states couphng to the quantum wire since bulk Landau states cannot be excited 
by the leads. The edge states shown in Figs. EJl tolHli correspond to the points in the transmission spectrum also 
labeled by (a) to (d) in Fig. [3 By comparison with the scattering wavefunctions (as in Fig. , we observe that in the 
magnetic field region B" < B < B"~^ edge states have up to n — 1 transverse nodes in the direction perpendicular to 
the boundary. Furthermore, the number of longitudinal antinodes from entrance to exit lead (see the corresponding 
numbers in Fig. |Hli,b) can be directly mapped onto successive maxima in Fig. [3 (see numbers there). The range of B 
depicted in Fig. corresponds to _B > B^ at kp = l.b'K/d. The transmission spectrum becomes increasingly complex 
as B is reduced or equivalently kp is increased (not shown). 

To determine the positions of the values i?" we consider the energy shift of Landau levels near the boundary. Bulk 
Landau levels are degenerate since their quantized energy = (n + 1/2)B /c is independent of their positition in 
space. This degeneracy is lifted if a Landau state is placed in the vicinity of the cavity wall: with decreasing distance 
to this boundary the energy of the state increases. Therefore the energies of edge states associated with the quantum 
number n lie above the asymptotic bulk value En. When the incoming electron is diffracted at the mouth of the 
entrance lead, only those edge states whose energy is below the Fermi energy can carry the flux. Due to the sharp 
edges at the junction between lead and quantum dot, all energetically accessible edge channels are populated. The 
blow-up of the scattering wavefunction near the lead mouth (Fig.^ highlights the diffractive edge scattering. This is 
a contrast to smooth edges where states in the lead could cross the lead junction adiabatically, i.e. without changing 
their state of quantization.'^— With sharp lead junctions however all edge states with quantum number n' < n will 
participate in transport up to a magnetic field where i?„ touches the Fermi energy. En — Ep, i.e. at the critical 
magnetic fields B^/c — Ep/{n+ 1/2). These threshold values are indicated by the dot-dashed vertical bars in Fig. l7l8l 
for B^ « 71.1, B^ « 118.4 and B^ « 355.3 for a lead width d = 0.25 and kp = l.bw/d. In our numerical data, both 
the position of these threshold values as well as their independence of the geometry are in excellent agreement with 
this prediction. The only exception is the critical magnetic field B^. Its value (355.3) lies slightly above the point 
where the transmission spectrum ceases (at B « 351.8) (see Fig. (S)). The reason for this deviation is the fact that 
the termination point of the spectrum is not determined by the magnetic field of the lowest bulk Landau level in the 
cavity (i.e. B^), but by the highest field at which the leads still carry flux. In the leads, however, the magnetic length 
does not satisfy the condition <C c? (at kp — l.dir/d : Ib ~ d/A.7). Contrary to the cavity, the lead wavefunctions 
still "feel" the constriction by the opposing wall. For this reason the threshold magnetic field values of the transverse 
lead states lie slightly below those of the bulk Landau levels, leading to a termination already below B'^. 

D. Multi-channel interferences 

The regular oscillations above Bl as well as the complex fluctuating pattern below B^ can be explained by a multi- 
channel scattering description. This model can be viewed as a generalizationSSi^ of a single-channel picture42i^ 
For this description to be applicable, the cavity of the dot has to have smooth boundaries and disorder must be 
absent. Under these circumstances the flux transported by edge states is conserved in the interior and changed only 
by diffractive scattering at lead junctions: At the junction, a fraction of the flux will exit through the lead while the 
remaining portion of the flux will continue to propagate along the boundary. The stationary scattering state can be 
viewed as the coherent superposition of repeated loops around the billiard. In order to translate this picture into an 
analytic expression we define amplitudes for transmission and reflection at the two lead junctions. We denote the 
amplitudes for transmission from transverse mode m in the entrance lead to the edge state in the dot with quantum 
number i by tmi. The amplitudes stand for transmission from edge state i in the dot to the transverse mode n in 
the exit lead. The amplitudes f^y(r-j) describe edge state reflection at the entrance (exit) lead from mode i to mode 
j. (The tilde signs serve to distinguish these amplitudes from the transport coefficients for the whole geometry tnm 
and r„m.) We further define the following matrices 

[^]^, = iy-e^"^^^"'''^^^'. [r']„j. (3.1) 

where Lj,Aj (LpA'j) denote the lengths L and areas A the edge state j covers from entrance to exit (from exit to 
entrance) of the dot. The areas A can be determined in gauge-invariant form, although the corresponding classical 
orbits are not necessarily periodici^ The transmission through the whole cavity tji — [T],y is then written as a 
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geometric series of matrices, 

T = f{l + R'R{l + R'R{l + ...)))f' 
= f If^iR'Ry] f 



\i=0 



= T{1 - R' R)-^T' . (3.3) 

Equation H3.3|) serves as a convenient starting point for the analysis of the transmission fluctuations. Consider first 
the regime B > B^, where only the lowest transverse edge state is excited. In this case Eq. (|3.3|l reduces to its scalar 
versioE>^fl*^ 

rptot _ I, |2 ^ KllPKllP fo 4^ 

' l-2Re[fiifne^'y] + |fii|2|fnP ' ^ ' 

with 7 ~ fci(Li + L[) — B{Ai + A[)/c. As expected, the fluctuations of \tii{B /c)\'^ are determined by an Aharonov- 
Bohm type phase 7. At fixed kp, the oscillation period is AB = 2Trc/Af'^. By Af'^ ^ Ai + A[ we denote the 
area which the edge state acquires in one revolution around the dot. Taking into account that the dynamically 
accessible area of the edge state is somewhat smaller than the geometric area, 

^tot ^ ^dot = 4 + TT (see Fig. 0, 

the prediction for the oscillation period is in excellent agreement with our numerical findings. Equation H3.4(l also 
explains why the oscillation period is increasing with increasing B (see Fig. IHJ. This explanation makes use of the 
somewhat counterintuitive fact, that for increasing magnetic field skipping orbits with fixed quantum number n have 
an increasing mean distance from the boundary.^^ Consequently, a larger B field implies a smaller enclosed area A\°^ 
and therefore a higher oscillation period AB. Furthermore Eq. (|3.4I) accounts for the fact that for most structures the 
successive maxima of T*°* reach unity42. (The small deviation from this rule of the stadium with 90° lead orientaion 
will be explained below). In addition to unitarity, + l^^nP — 1] we have for identical lead junctions Fn = r[^. 

(We call two junctions identical if the local environment of their lead mouths is the same and their respective distance 
is larger than a few wavelenths.) Above B^ scattering of an edge state at a junction is essentially a one-dimensional 
process, for which the probability for transmission from left to right has the same magnitude as vice versa. Identical 
lead junctions therefore also imply t'n = tn, provided that the two corners of the lead junction have the same shape. 
If and only if all of the three above conditions are fulfilled, Eq. (|3.4|l yields T*°' = 1 at the resonance condition 
7 = 27m, n g Z. Since for the two circle geometries and for the 180°-stadium the two lead junctions are identical, we 
indeed find in these cases that |tii(i?)p periodically reaches unity. On the other hand, when the leads are attached 
to the stadium at an angle of 90°, one lead is attached to the straight section while the other is attached to the 
semicircle. The local environment of the two lead mouths is in this case different (i.e. the lead junctions are not 
identical), for which reason our numerical results do not quite reach — 1, when the resonance condition is 

fulfilled for this geometry (see inset of Fig. [71i and Fig. |HJ). Another relation exists between the resonance condition 
and the behaviour of edge states. In the closed cavity an edge channel always encloses an integer number of flux 
quanta [BA/{(f)Qc) = m e N]. Therefore, the resonance condition is met whenever the energy of an edge state in the 
closed cavity crosses the Fermi edge4i 

One interesting feature of the transmission fluctuations in the single-channel regime of the circular dot {B > B\) 
is their invariance with respect to the lead orientation. The numerical results for the transmission probabilities of 
the circle with 180° and 90° lead orientation differ only at the tenth (!) decimal digit. This fact, as well as the 
observation, that in the case of the stadium billard the two lead orientations give different results, can again be 
explained by Eq. H3.4|l . The important point to note is that the interference phase (7 = fciL*°* — BA^°^ /c) does 
not change when changing the positions of the leads around the circle. Due to the rotational symmetry also the 
coefficients in, t'n, rn and f'^^ are the same for different angles between the leads. The same is thus true for the total 
transmission T*°* through the circular dot. The only exception to this rule occurs when the two leads are in close 
proximity to within a few wavelengths. In this case the transmission probability changes as compared to the results 
for the 180° and 90° lead orientation (not shown). 

The fluctuations in the regime B < B]. displayed in Fig.|7|can be analyzed with the help of a multi-channel scattering 
description. In the interval B^ 1^ B < B^ two channels corresponding to two edge states are open in the cavity and 
one channel in each of the leads. From the entrance to the exit lead mouth the two edge channels acquire the phases 
^ikiLi-iBAi/c ^ik2L2-iBA2/c j-espectively. Their interference at the exit lead will therefore give an oscillatory 
contribution to the total transmission T'°*(i?) = of the form cos^[B{Ai — A2)/{2c)]. For closer analysis we 

need to evaluate Eq. (|3.3|) which involves the inversion of 2 x 2 matrices. In the case of parallel lead orientation the 
corresponding expressions are simplified due to the fact that the phases acquired from entrance to exit lead and vice 
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versa are the same {An = and _L„ = L'J, 



ill = 



(3.5) 



with the abbreviated notation ipn — fcn^n — BAn/c. In Fig. EH we show one half-period of the beats in T*°*(_B) = 
for [tto < B{Ai — A2)/ (2c) < 7r(n + 1)], as calculated with Eq. (|3.5(l . The absolute square of the numerator 
(dashed line, A'') and denominator (dotted line D) of Eq. (|3.5|) display very similar oscillations, both in frequency 
and amplitude. However, since T*°* = N/D, a series of dips are superimposed on the term cos^[i?(Ai — A2)/{2c)] at 
the points where N and D have their common minima. To classify these dips (i.e. antiresonances) we make use of 
the fact)S^ that the unitarity of Eq. (|3.5|l can be satisfied by mapping the transport coefficients at the lead junctions 
(which are assumed to be identical) onto six independent parameters (r, e, (j), 1?, 0i, ^2), 



ill = ?ii = ^(l-r2)ee^[Wi+V')/2+'5], 
h2 = i'21 = y(]~72)(T37)e»[(02+^)/2+i^]^ 
hi^P'n = -[{l-e) + er] e'(f'^+'>\ 
fi2 = f'i2 = r2i - ^21 = (1 - r)^e{l-e) e*[('^i+^2)/2+^] 



(3.6) 



^^22 = -[e+ (1 - e)^] e 



Out of this set of parameters, only two (e,r) are physically relevant. The variable e represents the coupling of the 
incoming lead state to the edge state n = 1 in the interior and r e M is related to the reflection coefficient f'/j^ of an 
incident channel at the lead mouth, r'l^ — re^^'^'^^K Both quantities e, \r\ are restricted to the interval [0, 1]. With 
these terms the absolute square of tu [Eq. 1)3. 5|l ] can be considerably simplified;^ 



(l-r2)2 



sin^(77/2 + i?o + A) + r2 



2 ' 



(3.7) 



with <p = {<p2+ (j)2) - {fi + (t>i), r] = {if2 + h) + {'Pi + <t>i). r' = {\ - e) e-'-^/^ ^ ^ g#/2^ g ^ arg(r'), ^ -d + 5, a 
1 + r e^(''+2''), 13=1 + re-2*^ A = arg(^/a). The linewidth Eq is given by 



To 



2r'l3/a 

In the generic case of r 7^ 0, resonances occuring in Eq. (|3.7(l show a typical Fano profile of the forni^ 



rpt' 



bK|2 



{B/c~Bn/cf 



{B/c~Bn|c+^Y+Tl 



(3.8) 



(3.9) 



with being the coefficient for background scattering 4^ The Fano resonances at B/c = Bn/c — A will have an 
asymmetric lineshape unless A = (i.e. r = 0). This is however the case for the billiard systems we consider, since 
almost no refiection of incoming lead states takes place at the lead mouths, f'd k, 0, and therefore r « 0. Under this 
assumption Eq. (|3.7|l simplifies to 



sin^(77/2 + 7?o) 



(3.10) 



with linewidth Tq = (1— |r'p)/(2|r'|). This equation describes symmetric resonance lineshapes which can be identified 
as window resonances (also called Breit-Wigner dips/antiresonances) of the form 



{B/c-Bn/c)^ 
{B/c-Bn/cf + Tl 



(3.11) 



The physical picture resulting from this analysis is the following: In the magnetic field region Bl < B < Bl, where two 
edge states are present in the interior of the dot and one in each of the leads, the transmission probability shows large- 
scale oscillations intersected by sharp window resonances. The large-scale envelope function is given by 1/(1 -I- Eg). 
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Its maxima perfectly match with the roughly estimated term cos'^[B{Ai — A2)/{2c)] from above and can therefore 
be identified with the numbered points in Fig. each of which corresponds to an integer number of longitudinal 
antinodes in the wavefunction along the boundary (see Fig. . The antiresonances superimposed on these oscillations 
occur at magnetic fields B — Bn (where 77/2 + iJo = nir, n € 1) and their linewidth is given by Fq. As a result, 
resonances which are situated on maxima of the term 1/(1 + Fg) are sharper than at its mimima [see numerical data 
in Fig. 13 for confirmation]. For an increasing number of edge states populated in the cavity {B < B^) our numerical 
results show that the density of antiresonances is rapidly growing. This behaviour finally leads to a resonance overlap 
for a large number of edge states, which is prerequisite for the onset of Ericson fluctuations (i.e. universal conductance 
fluctuations). 

For completeness we remark that the above analysis for the B-dependence of T*°* can likewise be carried out with kp 
instead of B as the variable parameter. We can similarly identify threshold values fc" for kp, below which a number 
of n edge states survive. The numerical results for the transmission probablility T^°^{kp) for the case of one or two 
participating edge states (not shown) can again be described by Eq. 13.4(1 or Eq. (|3.7|l respectively. 



E. Comparison with experiments 

A series of experiment3i2iiii^2ii& have been performed where Aharonov-Bohm oscillations (ABOs) similar to the ones 
discussed here have been ob served in ballistic transport measurements. The origin of the ABOs in these experiments 
is however twofold: I n Refs . ll2l4lL it is the presence of edge states in a quantum dot which gives rise to the observed 
oscillations. In Refs. I47l48l on the contrary, the investigated scattering devices have the form of a ring, to which the 
scattering wave function is confined. The latter setup thus gives rise to ABOs already at low magnetic fields and 
has therefore been more readily accessible to a theoretical description!^ However, to our knowledge, no quantitative 
description for magnetotransport through a quantum dot in the regime of only one or two participating edge states 
has yet become available. We therefore discuss in the following similarities and differences between the experimental 
data and our calculations in this field. One important observation is that the magnetic fields where these quasi-regular 
transmission fluctuations appear in the experiment are lower than in the present calculation. For example, in the 
experiment for circle and stadium shaped quantum dots in a Ga^s/^/Ga^s-heterostructureji^ the threshold magnetic 
field values would be (in Sl-units) 

27r/i , , /27r 

Bc = to ^iu2 ^ith Xp^J—. 3.12 

With a given sheet density of Ug = 3.6 x lO^^cm"^ in the interior of the dot, the threshold magnetic fields are given 
by B^ « 3 Tesla and B^ « 5 Tesla. However, in the experiment regular oscillations were already observed below 2 
Tesla. At those field values we find highly irregular transmission fluctuations corresponding to a threshold magnetic 
field i?" with n ^ 1, indicative of a high density of resonances and Ericson fiuctuations. We expect the origin of this 
discrepancy to lie in the absence of sharp edges in the experiment and, hence, of diffractive edge scattering. In the 
experimental quantum dot, the edges should be fairly smooth, leading to near-adiabatic transitions to edge states at 
the entrance to the quantum dot. Therefore fewer edge channels are excited than by diffractive edge scattering, where 
all energetically accessible channels up to n are populated. Our present results suggest that the observed transmission 
fiuctuations are a direct measure of the sharpness of the edges at the lead mouth. Therefore, investigations of quantum 
dots with varying sharpness of edges would be desirable. Since these are, however, difficult to fabricate we point to 
a different experimental approach, which is based on the analogy between transport in the edge state regime and 
field-free transport through a rectangle where only few propagating modes participate. Such structures are accessible 
for microwave experiments. ^iSSi The measured transmission through such a microwave device could provide a stringent 
test for the multi-channel interference model presented above. 



IV. SUMMARY AND OUTLOOK 



We have presented a new technique for calculating ballistic magnetotransport through open quantum dots. The 
Modular Recursive Green's Function Method (MRGM) is an extension of the widely used standard recursive Green's 
function technique and is based on the decomposition of non-separable scattering geometries into separable substruc- 
tures (modules) . An unprecedented energy and magnetic field range can thereby be explored with high accuracy. We 
applied the MRGM to transport coefficients and scattering wavefunctions in the two extreme cases of high magnetic 
fields and short wavelengths. For very small cyclotron radii we found periodic oscillations in the transmission spectrum 
and beating phenomena, which are restricted to well defined intervals (as a function of B and kp likewise). These 
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features could be explained by interferences between edge states, travelling along the boundary of the cavity. For 
these states scattering only takes place at the lead junctions, whose sharp edges play a crucial role for the dynamics 
of the system. For a detailed analysis a multi-channel interference model was employed. This model allows to classify 
the observed transmission fluctuations in the framework of Fano resonances. For only one edge state present in the 
circular dot transport is independent of the lead orientation provided that the lead mouths are identical and separated 
from each other. Future envisioned applications include the investigation of Andrccv billiards^ quantum Hamilto- 
nian ratchets, fractal conductance fluctuations, ^■^•^'^ and shot noise. Furthermore, the MRGM also seems suitable 
to perform ab-initio calculations of the integer Quantum Hall effect. '"^^ The challenge is in this case the inclusion of a 
disorder potential which is compatible with the separability conditions. 
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FIG. 1: Joining and disconnecting of modules by application of a Dyson equation: two semi-infinite leads. The hard wall 
boundary conditions at the sites on the border of the modules are represented by empty circles (accessible space by full circles). 
The gray shaded areas P and Q are those grid slices at which the Green's functions are evaluated (see text). 




FIG. 2: Applying a Dyson equation in "reversed mode" to construct Green's functions for (a) a rectangle out of a cylinder 
surface and (b) a semi-circle out of a full circle, respectively. In (a) the periodic boundary conditions are transformed into hard 
wall boundary conditions. The gray shaded areas P, Q and X are those grid slices at which the Green's functions are evaluated 
(see text). 
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FIG. 3: Applying a Dyson equation to construct Green's functions for (a) a semi-infinite lead and (b) a stadium billiard out 
of "modules". In (a) joining a transverse slice with a semi-infinite lead schematically leaves the Green's function of the lead 
invariant. In (b) an additional link module is added to facilitate the coupling between the half-circle and the rectangle module. 
Notation as in Figs. and |5| 
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FIG. 4: Comparison between the present MRGM (solid line) and the wavefunction matching technique^ (dotted line) for the 
first-mode transmission probability |tii(fcF)|'^ at B/c = 1 in a small window of kp- (a) circle with perpendicular leads, (b) 
stadium with perpendicular leads (d = 0.35, A*^"* =4-1- tt). In both cases also |tii(fcF)P -I- |rii(fcF)P is shown. Contrary to the 
MRGM (solid line), the wave function matching technique (dotted line) deviates from the unitarity limit in (b). 
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FIG. 5: (color online) Absolute square of the scattering wave functions at high kp [(a): kp = 257r/d, (b): kp = 

12.57r/ci, (c),(d): kp ~ 6.0l7r/ci] for the four quantum dots considered: circle and stadium with relative lead orientation of 90° 
and 180°, area ^4^^°' = 4 + tt and lead width d — 0.25. The localization around classical trajectories (see insets for comparison) 
is clearly visible. In Figs, (a)-(c) the magnetic field B — 0. In Fig. (d) the magnetic field B/c — 30.5 allows for a whole bundle 
of equivalent trajectories with cyclotron radius Tc = kpc/B ~ 2.48 to contribute to transport. 



FIG. 6: (color online) Absolute square of the scattering wave functions \4>{x,y)f in the edge state regime. The area of all 
geometries A"^"^ = 4 + tt, lead width d = 0.25, and kp = 1.5'K/d. The four plots correspond to the points in the transmission 
spectra (Fig. 0, indicated by the letters (a)-(d). The numbers along the longitudinal direction of the edge states count the 
number of antinodes between entrance and exit lead (see corresponding numbers in Fig. |7|l. Note that edge states in the 
magnetic field region B"^^ < B < B" have up to n transverse nodes: (a) circle, 180°, n — 2, (b) circle, 90°, n = 2, (c) stadium, 
180°, n = 3, and (d) stadium, 90°, n — 1. 
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FIG. 7: Transmission probabilities \ti\{B/c)'[^ in the high magnetic field limit for circle and stadium billiard with 180° or 
90° lead orientation, [kp ~ l.!3Tv/d,d = 0.25). Bl and are the threshold magnetic fields B"/c — kp/{2n + 1) (vertical 
dash-dotted lines). Above BI regular oscillations appear (see insets for magnification). For B^ < B < B^ irregular fiuctuations 
set in. Their large-scale structure can be explained by the number of interference maxima the two edge states form along the 
boundary between entrance and exit lead (see indicated numbers). The points (a)-(d) correspond to the wavefunctions shown 
in Fig. il 




FIG. 8: Transmission probabilities |iii(B/c)p in the high-field limit, near the point where transport terminates. The dotted 
line stands for the circle billiards in both lead geometries (their transmission probabilities axe identical) and the dashed line 
for the stadium billiard with 90° lead geometry. The solid curves represent the upper and lower bounds of the oscillations 
(offset for better visibility) . The two dash-dotted vertical lines mark the point where transport breaks off and the analytically 

determined threshold value ~ 355.3 (see text for details). The inset shows that the transmission probabilities for the circle 
reach the maximum value 1 which is only approximately true for the stadium with 90°-lead geometry. 



FIG. 9: (color online) Electron density \'ijj{x,y)\^ for the circle billiard with diffractive scattering highlighted. {A'^°^ — 4 + n, 
lead width d = 0.25 and kp = 1.57r/d.) The magnetic field B = 118.7 is just above the threshold to the single-edge state regime 
Bl = 118.44. 
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interference phase 

FIG. 10: One half-period of the beating nn < B{A\ — A2)/{2c) < 7r(n + 1) in the transmission probabihty |tii(_B/c)p (sohd 
hne), as calculated with our interference model [see Eq. 13.7t 1. The nominator (dashed line, A'^) and denominator (dotted line, 
D) of |tii(B/c)p show very similar oscillations (with a small offset). (N/D) features sharp "dips", at the points where and 
D have their common minima. These dips are window resonances (also called Breit-Wigner antiresonarices) and represent a 
symmetric limit of the Fano resonance lineshape. See text for details. 
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